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Abstract 

The effects of decoherence for quantum system coupled with a bosonic field are investi- 
gated. An application of the stochastic golden rule shows that in the stochastic limit the 
dynamics of such a system is described by a quantum stochastic differential equation. 

The corresponding master equation describes convergence of a system to equilibrium. 
In particular it predicts exponential damping for off-diagonal matrix elements of the sys- 
tem density matrix, moreover these elements for a generic system will decay at least as 
exp(— tN^), where N is a number of particles in the system. 

As an application of the described technique a derivation from first principles (i.e. starting 
from a Hamiltonian description) of a quantum extension of the Glauber dynamics for systems 
of spins is given. 

1 Introduction 

In the present paper we investigate a general model of quantum system interacting with a bosonic 
reservoir via an Hamiltonian of the form 



H = H + \Hj 



where Hq is called the free Hamiltonian and Hj the interaction Hamiltonian. 

The stochastic golden rules, which arise in the stochastic limit of quantum theory as natural 
generalizations of Fermi golden rule [|TJ, 0, provide a natural tool to associate a stochastic flow, 
driven by a white noise equation (stochastic Schrodinger) equation, to any discrete system inter- 
acting with a quantum field. This white noise Hamiltonian equation which, when put in normal 
order becomes equivalent to a quantum stochastic differential equation. The Langevin (stochas- 
tic Heisenberg) and master equations are deduced from this white noise equation by means of 
standard procedures which are described in Q. 

We use these equations to investigate the decoherence in quantum systems. 

In the work || , extending previous results obtained with perturbative techniques by , it was 
shown on the example of the spin-boson Hamiltonian that the decoherence in quantum systems 
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can be controlled by the following constants (cf. section 2 for the definition of the quantities 
involved) 



For the simplest case of the equalibrium state of the reservoir with the temperature = kT this 
constant will be equal to (actually this is true for large temperatures and for the dispersion 
function u(k) = \k\). 

In this paper we extend the approach of from 2-level systems to arbitrary quantum systems 
with discrete spectrum. Our resuts show that the stocastic limit technique gives us an effective 
method to control quantum decoherence. 

We find that, under the above mentioned interaction, all the off-diagonal matrix elements, of 
the density matrix of a generic discrete quantum system, will decay exponentially if Ke(g\g) are 
nonzero. In other words we obtain the asymptotic diagonalization of the density matrix. 

Moreover, we show that for generic quantum system the off-diagonal elements of the density 
matrix decay exponentially as exp(— iVRe (g\g)t), with the exponent proportional to the number 
iV of particles in the system. Therefore for generic macroscopic (large N) systems the quantum 
state will collapse into the classical state very quickly. This effect was built in by hands in several 
phenomenological models of the quantum measurement process. In the stochastic limit approach 
it is deduced from the Hamiltonian model. 

This observation contributes to the clarification of one of the old problems of quantum theory: 
Why macroscopic systems usually behave classically? i.e. why do we observe classical states 
although the evolution of the system is a unitary operator described by the Shrodinger equation? 

Moreover, this result allows to distinguish between macroscopic systems (that behave classi- 
cally) and microscopic systems (where quantum effects are important). Quantum effects (or ef- 
fects of quantum interference) are connected with the off-diagonal elements of the density matrix. 
Therefore the following notion is natural: the macroscopic system is a system where off-diagonal 
elements of the density matrix decay quickly (faster than the minimal time of observation) . Us- 
ing that off-diagonal elements decay as exp(— NRe (g\g)t), we get the following definition of the 
macroscopic system: NRe(g\g) » 1. 

The quantum Markov semigroup we obtain lives invariant the algebra of the spectral pro- 
jections of the system Hamiltonian and the associated master equation, when restricted to the 
diagonal part of the density matrix, takes the form of a standard classical kinetic equation, describ- 
ing the convergence to equilibrium (Gibbs state) of the system, coupled with the given reservoir 
(quantum field). 

Summing up: the convergence to equilibrium is a result of quantum decoherence. 

If we can control the interaction so that some of the constants Re(g\g) are zero, then the 
corresponding matrix elements will not decay in the stochastic approximation, i.e. in a time scale 
which is extremely long with respect to the slow clock of the discrete system. In this sense the 
stochastic limit approach provides a method for controlling quantum coherence. 

The general idea of the stochastic limit (see JTJ) is to make the time rescaling t — > t/X 2 in 
the solution of the Schrodinger (or Heisenberg) equation in interaction picture = e ttH o e ~ itH ; 
associated to the Hamiltonian H, i.e. 
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with Hj(t) = e ltH °Hie ttH °. This gives the rescaled equation 




(1) 



and one wants to study the limits, in a topology to be specified, 





(2) 



The limit A — » after the rescaling t — > t/A 2 is equivalent to the simultaneous limit A — > 0, t — > oo 
under the condition that A 2 t tends to a constant (interpreted as a new slow time scale). This 
limit captures the dominating contributions to the dynamics, in a regime of long times and small 
coupling, arising from the cumulative effects, on a large time scale, of small interactions (A — > 0). 
The physical idea is that, looked from the slow time scale of the atom, the field looks like a very 
chaotic object: a quantum white noise, i.e. a 5-correlated (in time) quantum field b*(t,k), b(t,k) 
also called a master field. 

The structure of the present paper is as follows. 

In section 2 we introduce the model and consider its stochastic limit. 

In section 3 we derive the Langevin equation. 

In section 4 we derive the master equation for the density matrix and show that for non-zero 
decoherence the master equation describes the collapse of the density matrix to the classical Gibbs 
distribution and discuss the connection of this fact with the procedure of quantum measurement. 

In section 5, using the characterization of quantum decoherence obtained in section 4 and 
generalizing arguments of , we find that our general model exhibits macroscopic quantum effects 
(in particular, conservation of quantum coherence). These effects are controllable by the state of 
the reservoir (that can be controlled by filtering). 

In section 6 we apply our general scheme to the model of a quantum system of spins interacting 
with bosonic field and derive a quantum extension of the Glauber dynamics. 

Thus our stochastic limit approach provides a microscopic interpretation, in terms of funda- 
mental Hamiltonian models, to the dynamics of quantum spin systems. Moreover we deduce the 
full stochastic equation and not only the master equation. This is new even in the case of classical 
spin systems. 

2 The model and it's stochastic limit 

In the present paper we consider a general model, describing the interaction of a system S with 
a reservoir, represented by a bosonic quantum field. Particular cases of this general model were 
investigated in ||, 0, |J. The total Hamiltonian is 



acting in the representation space T corresponding to the state (•) of bosonic reservoir generated 
by the density matrix N that we take in the algebra of spectral projections of the reservoir 



H = H + XH I = H S + H R + XH! 



where Hr is the free Hamiltonian of a bosonic reservoir R: 
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Hamiltonian. The reference state (•) of the field is a mean zero gauge invariant Gaussian state, 
characterized by the second order correlation function equal to 

(a{k)a*{k')) = (N(k) + l)5(k - k') 

(a*{k)a{k')) = N{k)5{k - k') 

where the function N(k) describes the density of bosons with frequency k. One of the examples 
is the (gaussian) bosonic equilibrium state at temperature 

The system Hamiltonian has the following spectral decomposition 

H s = Y j e r Pe T 

r 

where the index r labels the spectral projections of Hs- For example, for a non-degenerate 
eigenvalue e r of Hs the corresponding spectral projection is 

P £r = \e r )(e r \ 

where \e r ) is the corresponding eigenvector. 

The interaction Hamiltonian Hj (acting in Tis <8> J~) has the form 

Hj = E (D* ® A( gj ) + Dj ® A*( 9j )) , A(g) = J dkg(k)a(k) , 

3 

where A(g) is a smeared quantum field with cutoff function (form factor) g(k). To perform the 
construction of the stochastic limit one needs to calculate the free evolution of the interaction 
Hamiltonian: H^t) = e ftH °fl/e -itfl ' '. 
Using the identity 

r 

we write the interaction Hamiltonian in the form 

Hi = E E Pe„D*P £ , J dkg 3 (k)a(k) + h.c. (3) 

j rr' 

Let us introduce the set of energy differences (Bohr frequencies) 

F = {uj = e r — e r i : e r , e r > G Spec H s } 
and the set of all energies of the form 

F u = {s r '■ 3e r ' (e r , e r i e Spec Hs) such that e r — e r > = u} 
With these notations we rewrite the interaction Hamiltonian (|3|) in the form 

ff/ = EE E Pe T D*P £r ^ f dk5j(k)a(k) + h.c. = 

= Y,HK ( D j) I dkg 3 (k)a(k) + h.c. (4) 

j u£F 



4 



lim - exp - f(k) ) = 2*6(t)6(f(k)) (7) 



where 

E U (X) := £ P £r -.XP £r (5) 

It is easy to see that the free volution of E U (X) is 

e^ s ^(X)e-^ s = e-^E^X) 
Using the formula for the free evolution of bosonic fields 

e itHR a(k)e- itHR = e~ ituj(k) a(k) 
we get for the free evolution of the interaction Hamiltonian: 

= E E K {Dj) I dkg j (k)e* u W- u )a(k) + h.c. (6) 
j caeF 

In the stochastic limit the field Hi(t) gives rise to a family of quantum white noises, or master 
fields. To investigate these noises, let us suppose the following: 

1) u(k) > 0, Vfc; 

2) The d — 1-dimensional Lebesgue measure of the surface {k : us(k) = 0} is equal to zero (so 
that 5(u(k)) = 0) (for example u(k) = k 2 + m with m > 0). 

Now let us investigate the limit of Hj(t/X 2 ) using one of the basic formulae of the stochastic 
limit: 

1 fit 

^ exp b 

which shows that the term S(f(k)) in (0) is not identically equal to zero only if f(k) = for some 
k in a set of nonzero d — 1-dimensional Lebesgue measure. This explains condition (2) above. 

The rescaled interaction Hamiltonian is expressed in terms of the rescaled creation and anni- 
hilation operators 

ax,u,(t, k) = \ e~*& M^aik), lu £ F 
A 

After the stochastic limit every rescaled annihilation operator corresponding to any transition 
from e r i to e r with the frequency u = e r — e r > generates one non-trivial quantum white noise 

b u (t, k) = lim a^(t, k) = lim \ ^ a(Jfe) 

with the relations 

[b u (fi,k)XV,V)] = \™[axAt,k),a* x Jt',k')] = 

A — >U 

= lim e- iS fr {w{k) -^5{k - k') = 2n5(t - t')8(u(k) - u)5(k - k') (8) 

[bUt,k),b* ujl (t',k')] = 

(cf. ^). This shows, in particular that quantum white noises, corresponding to different Bohr 
frequencies, are mutually independent. 

The stochastic limit of the interaction Hamiltonian is therefore equal to 

Kt) =EEC ( D i) f dkg^b^t, k) + h.c. (9) 

j LUtF 
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The state of the master field (white noise) b w (t,k), corresponding to our choice of the initial 
state of the field, is the mean zero gauge invariant Gaussian state with correlations: 

(6* (t, k)bu(t', k')) = 2Tx5{t - t')8{u(k) - u)5{k - k')N(k) 

(b u (t, k)bl{t', k')) = 2n5(t - t')5{uj{k) - u)5{k - k')(N(k) + 1) 

and vanishes for noises corresponding to different Bohr frequences. 

Now let us investigate the evolution equation in interaction picture for our model. According 
to the general scheme of the stochastic limit, we get the (singular) white noise equation 

jU t = -ih(t)U t (10) 
whose normally ordered form is the quantum stochastic differential equation 

dU t = (-idH(t) - Gdt)U t (11) 

where h(t) is the white noise Hamiltonian @ given by the stochastic limit of the interaction 
Hamiltonian and 

dH{t) = E E (K Pi) dB JU (t) + E u (Dj) dB^(t)) (12) 
j u>eF 

r ft+dt 

dB iu {t) = j dkg^k) j 6 w (r,fc)dr (13) 

According to the stochastic golden rule (|TT|) the limit dynamical equation is obtained as follows: 
the first term in ( |TT| ) is just the limit of the iterated series solution for (|I|) 

l r t+dt ( r \ 

The second term Gdt, called the drift, is equal to the limit of the expectation value in the reservoir 
state of the second term in the iterated series solution for ([]]) 



if^/w^HI)) 



lim — 

A^O A 

Making in this formula the change of variables r = t 2 — ti we get 

1 ft+dt rO /f \ fl 



lim , 

A^O A 2 



Computing the expectation value and using the fact that the limits of oscillating factors of the 
form lim vanish unless the constant c is equal to zero, we see that we can have non-zero limit 

icti 

only when all oscillating factors of a kind (with ti) in (fTJ) cancel. In conclusion we get 



G = E E [° dr ([ dkg t {k) gj {ky<^-"\N{k) + l)E* (A) E u (Dj) + 

ij u>£F J -°° KJ 

+ f dk gi (k)^e^ T M k WN{k)E u (Di)E*(Di 
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and therefore, from the formula 



e^dt = % — = 7i5(u) - % P.P. - (15) 

-oo UJ — ZV UJ 



we get the following expression for the drift G: 



= EE (tela-CK (a) (£),•) + (A) c (A)) («) 

Let us note that for ( 16|) we have the following Cheshire Cat effect found in ||: even if the 
frequency u is negative and therefore does not generate a quantum white noise the corresponding 
values {g\g)u in (|16|) will be non-zero. In other terms: negative Bohr frequencies contribute to an 
energy shift in the system, but not to its damping. 

Remark If F is any subset of Specifj and X r are arbitrary bounded operators on Hs then 
for any t G R 

e uH s P £r X r P £r = £ e u ^P £r X r P £r = Pe r X r P £r e itHs 
e r eF e r eF e r eF 

In other words: ^ £t&f P Er X r P er belongs to the commutant L°° (Hg) of the abelian algebra 
L°° (Hs), generated by the spectral projections of Hs- 

A corollary of this remark is that, for each uj G F, for any bounded operator X G L°° (Hs) 
and for each pair of indices (i,j) the operators 

E u (Di)XE*(Dj), K(Di)XE u (Dj) (17) 

belong to the commutant L°° (Hs)' of L°° (Hs)- In particular, if Hs has non-degenerate spectrum 
so that L°° (Hs) is a maximal abelian subalgebra of B(Hs), the operators (|i~7j) also belong to 
L°°(H S ). 



3 The Langevin equation 

Now we will find the Langevin equation, which is the limit of the Heisenberg evolution, in inter- 
action representation. Let X be an observable. The Langevin equation is the equation satisfied 
by the stochastic flow j t , defined by: 

j t (x) = u;xu t 

where U t satisfies equation ([□]) in the previous section, i.e. 

dU t = (-idH(t) - Gdt)U t (18) 
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To derive the Langevin equation we consider 

dj t (x) = j t+ dt{x) - j t (x) = du;xu t + u*xdu t + du;xdu t (19) 

The only nonvanishing products in the quantum stochastic differentials are 

dBUt)dB* 3UJ {t) = 2Re (g t \g 3 )Zdt, dB* w (t)dB JUJ (t) = 2Re (g t \g^dt (20) 
Combining the terms in (|19|) and using flT8|), ([T2]), ([16]) and ( p0| ) we get the Langevin equation 

dj t (x) = Y,jt o e a (x)dM a (t) = Yl Jt ° e nju ,(x)dM nj "(t)+j t o e (x)dt (21) 

where 

dM- l ^{t) = dB ju (t), 0-ij U (X) = -i[X, El (Dj)] (22) 
dM^(t) = dB* ju> (t), O^uiX) = -i[X, A (Dj)] (23) 

and 

0o(X) =EE(- J Im K (A) A (A)] + * Im (ftbi)^^, A (A) K (A)]+ 

+2Re (g t \g 3 )Z (k (A) *A (Dj) - ~{X, E* (A) A (A)}) + 

+2Re (e u (A) x^; (A) - \{X, A (A) A (A)})) (24) 

is a quantum Markovian generator. The structure map 6 (X) has the standard form of the 
generator of a master equation |7|] 

9 (X) = vl/ (X) - ±{*(l),X} + i[H,X] 

where \l/ is a completely positive map and H is selfadjoint. In our case ^(X) is a linear combination 
of terms of the type 

C(A)XA(A) 

Remark A corollary of the remark at the end of section 2 is that the Markovian generator # 
maps L°°(Hs)' into itself. Moreover, if X in (|2~3| ) belongs to the L°°(Hs) then the Hamiltonian part 
of 6q(X) vanishes and only the dissipative part remains. In particular, if Hs has non-degenerate 
spectrum then 8q(X) maps L°°(Hs) and has the form 

0o(X) =EE ( 2 Re (9i\9j)Z (A (A) XE U (Dj) - XEl (A) A (A)) + 

+2Re ( 9i \ gj )+ (A (A) XEl (A) - XE U (A) A (A))) 
for any X G L°°(H S ). 

The structure maps 6^ in (|2ll) satisfy the following stochastic Leibnitz rule, see the paper || . 
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Theorem. For any pair of operators in the system algebra X , Y, the structure maps in the 
Langevin equation ( pip satisfy the equation 



e a {XY) = e a {x)Y + xe a {Y) + Y,c^e p {x)e 1 {Y) 

where the structure constants c^ 7 is given by the Ito table 

dM p {t)dM\t) = Y,<%dM a {t) 



The conjugation rules of dM a (t) and 9 a are connected in such a way that formula fl2J|) defines a 
*-flow (* o jt = jt° *)■ 

3.1 Evolution for the density matrix 

Let us now investigate the master equation for the density matrix p. 

We will show that if the reservoir is in the equilibrium state at temperature f3~ l then for the 
generic system with decoherence the solution of the master equation p(t) with t —>■ oo tends to 
the classical Gibbs state with the same temperature f3~ l . This phenomenon realizes the quantum 
measurement procedure — the quantum state (density matrix) collapses into the classical state. 

To show this we use the control of quantum decoherence that was found in the stochastic 
approximation of quantum theory, see and discussion below. 

Let us consider the evolution of the state (positive normed linear functional on system ob- 
servables) given by the density matrix p, p(X) = tr pX. The evolution of the state is defined as 
follows 

Pt = ftip) = P°3t 
Therefore from (pT|) we get the evolution equation 



dpt{X) = p o d 3t (X) = po^j t o 9 a (X)dM a (t) = Y,Pt (9 a (X)dM a (t)) 

a a 

Only the stochastic differential dt in this formula will survive and we get the master equation 

j t p t (x) = p t oe (x) = e;(p t )(x) (25) 

Let us consider the density matrix p = ps <S> Pr, 

Ps,t = J2p(p, v,t)\fi)(v\ 

where \p), \v) are eigenvectors of the system Hamiltonian Hg- 
Using the form ([24|) of 9 and the identities 

tr Y[X,A] = -tr %A]X 
tr Y (AXB - ~{X, AB}\ = tr (bYA - -{Y, AB}j X 
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the master equation (^Sj) will take the form 

E T+ p ^> ^ 1^ ^1 = E ^ *) 
EE V * Im (ftb)J[l^)H^:(A)^(^)] -< Im (g i \g J )Z[\iJ.)(v\,E u (D i )K(Dj)]+ 

ij wgF ^ 

+2Re( gi \ 9j )Z (e u (Dj) \fj)(u\EZ{D t ) - ^{\n)(u\,E^(Di)E u + 
+2Re (iC (Dj) |//)(^|^ (A) - (A) K (A)})) = 

= E Im fa\9i)Z (\l*)(v\x u {e»)DiP<*-+,D j P €y - + 

+2Re(^|^)j^M^(^)^ M -^jlM)(^lA* P ^-^- 

-i Im (gi\gj)+ [\n)(u\x-u(e v + u)DiP £i/+u) D*P £u - X-u>{e^ + u)P £fl DiP £fi+LU D*\fi)(u\j + 

+2Re (ftbj)w (x-w(e/i + ^)X-^ + w)P e(t+w D^|^)<i/|A-P ev +«- 

(|a*X^IX-^(^ + uj)D % P £v+u1 D*P £u + x-.^ + W )P e(1 A^ + ^|M}(H))) (26) 
where x^(^m) = 1 if ^ -^L an d equals to otherwise. 

4 Dynamics for generic systems 

Let us investigate the behavior of a system with dynamics defined by fl26|) . This dynamics will 
depend on the Hamiltonian of the system. 

We will call the Hamiltonian Hg generic, if: 

1) The spectrum Spec Hs of the Hamiltonian is non degenerate. 

2) For any Bohr frequency to there exists a unique pair of energy levels e, e' G Spec Hs such 
that: 

uj = e — e' 

We investigate (|26| ) for generic Hamiltonian. We also consider the case of one test function 
gi(k) = g(k), although this is not important. In this case 

E U (X) = \a')(a'\X\a)(a\ = \a')(a\(a'\X\a) 

where u = e a — e a > . The Markovian generator 9q in (p6|) takes the form 

a,a> V 
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+2Re{g\g)- al {\a'){a'\{a\X\a)-\{X^a){a\}y 

-i Im (£b)+ CT ,[Xj<^ (27) 
We use here the notion 

Notice that the factors Re (g\g) aa i are > only for e a > e a / and vanish for the opposite case. 
It is easy to see that the terms in fl27|) of the form 

\a)(a\(a'\X\a') 

for off-diagonal elements of the density matrix X = \p)(v\ are equal to zero. We will show that 
in such case the equation fl26|) will predict fast damping of the states of the kind 

In the non-generic case one can expect the fast damping of the state with different 

energies and e u . 

With the given assumptions the action of 9q on the off-diagonal matrix unit |/i)(z/|, ^ e v is 
equal to A^p) where the number A^ v is given by the following 

+i\m{g\g)UW\a)\ 2 -Re{g\g)^\{a\DW 
-Re(0b)-|<*^ (28) 

The map 6q multiplies off-diagonal matrix elements of the density matrix p$ by a number A^ u . 
Let us note that 

ReV <0 

Moreover, for generic Hamiltonian the map 9$ mixes diagonal elements of the density matrix 
but does not mix diagonal and off-diagonal elements (the action of 8$ on diagonal element is equal 
to the linear combination of diagonal elements). 

The equation (^6|) for the generic case takes the form 

+ £ E(pK^) (2Re {g\g)-, a \(o\D\o>)\* + 2Re (^1^)+^ | (^'|^k> | 2 ) - 

-p(a,t)(2Re(^|^|(a|D|aO| 2 + 2Re(^|^|(a'| J D|a)| 2 )) (29) 

with A^y given by ( P8"D and p(cr,t) = p(a,cr,t). 
For instance we get 

= exp{A^t)\p)(v\ 

We see that if any of Re (f7|(7)|(/3|ZP|c^)| 2 in (|28|) is non-zero then the corresponding off-diagonal 
matrix element of the density matrix decays. We obtain an effect of the diagonalization of the 



11 



density matrix. This gives an effective criterium for quantum decoherence in the stochastic ap- 
proximation: the system will exhibit decoherence if the constants Re {g\g^ are non-zero. 

Now we estimate the velocity of decay of the density matrix for a quantum system with 

TV particles. The eigenstate \p) of the Hamiltonian of such a system can be considered as a tensor 
product over degrees of freedom of the system of some substates. Let us estimate from below the 
number of degrees of freedom of the system by the number of particles that belong to the system 
(for each particle we have few degrees of freedom). To get the estimate from below for the velocity 
of decay we assume that |(a|.D|/i)| 2 in fl2"S| ) is non-zero only if the state a differs from the state p, 
only for one degree of freedom. 

Then the summation over uj (or equivalently over a) in fl28|) can be estimated by the summation 
over the degrees of freedom, or over particles belonging to the system. If we have total decoherence, 
i.e. all Re (g\g) are non-zero, then, taking all corresponding |(cr|Z)|/x)| 2 = 1, we can estimate ( p8f) 
as — iVRe (g\g), where N is the number of particles in the system, or 

= exp(-iVRe(<7|<7)0|/i><H (30) 

The off-diagonal element of the density matrix decays exponentially, with the exponent propor- 
tional to the number of particles in the system. Therefore for macroscopic (large N) systems with 
decoherence the quantum state will collapse into the classical state very quickly. 

This observation clarifies, why macroscopic quantum systems usually behave classically. The 
equation (|30| ) describes such type of behavior, predicting that the quantum state damps at least 
as quickly as exp(— iVRe (g\g)t). Therefore a macroscopic system (large N) will become classical 
in a time of order (iVRe {g\g))^ 1 . 

Let us estimate the constant Re (g\g) for the equilibrium state of the reservoir with the tem- 
perature = kT. In this case 

Re (9\9)Z = l^-r-, I dk \g(k)\ 2 5(uj(k) - U ) 



h 1 - e-#" 

Taking g{k) = 1 and using that the dispersion function u(k) depends only on \k\ we get 



dk\g(k)\ 5(u(k)-u)= / dVt I dpd(uj(p) - u) = Air / dpS(u(p) - to) 

J Jo Jo 

where / dfl is the integration over angles. If we take the dispersion function ui{k) = \k\, then for 
this integral we get Akuj. 

Therefore for this choice of dispersion function we get 

Re(g|o),~T = — — 

and analogously 

/ i \ _l 1 Att 2 uj 
Re (g o),, = ——3 

In the limit of small (3 (or high temperatures) f3~ x = kT » uj both these integrals tend to 

A 2 kT 

An — 
n 
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Summing up, we get that the constant Re (g\g)^ for the case of high temperature will be equal 
to up to multiplication by a dimensionless constant depending on the model. 

This means that every degree of freeedom that energetically admissible (kT » uo) and not 
forbidden by the model (|(^|-D|/i)| 2 7^ 0) gives the term of order ^ in the exponent for dumping 
of off-diagonal matrix elements. 

The off-diagonal matrix element will dump as exp(— tN^), where N is the number of degrees 
of freedom (that for a generic system can be taken proportional to the number of particles). Off- 
diagonal matrix elements describe the quantum interference. Our result for the dumping of off- 
diagonal matrix elements (|30D gives us a possibility to distinguish between microscopic systems 
(where quantum effects such as quantum interference are important) and macroscopic system 
which can be described by classical mechanics. The macroscopic system is a system satisfying 

kT 

N— » 1 (31) 

Actually the value N^- is of dimension of t~ l } and fl3~ij ) means that (N^jpj is much less than 
the time of observation. 

In the last section of the present paper we will illustrate the collapse phenomenon ( ftOD using 
the quantum extension of the Glauber dynamics for a system of spins. 

We see that the stochastic limit predicts the collapse of a quantum state into a classical state 
and, moreover, allows us to estimate the velocity of the collapse (j30|) . One can consider ( p0|) as 
a more detailed formulation of the Fermi golden rule: the Fermi golden rule predicts exponential 
decay of quantum states; formula ( j50| ) also relates the speed of the decay to the dimensions 
(number of particles) of the sstem. 

Consider now the system density matrix p$ G C, where C is the algebra generated by the 
spectral projections of the system Hamiltonian Hs, and consider the master equation ( p7|) (we 



consider the generic case). We will find that the evolution defined by this master equation will 
conserve the algebra C and therefore will be a classical evolution. We will show that this classical 
evolution in fact describes quantum phenomena. 

For ps,t G C we define the evolved density matrix of the system 

Ps,t = J2p(cr,t)\a)(a\ 

a 

For this density matrix the master equation (f25|) takes the form 



j t p(a,t) (2Re (g\g)- l(7 \(a\D\a')\" + 2Re (g\g)i a ,\{*'\D\a)\ 



-p(a,t)(2Re(^|^ ff |(a|D|a')| 2 + 2Re(^|^|(a'| J D|a)| 2 )j (32) 
Let us note that if p(a, t) satisfies the detailed balance condition 

p(a,t)2Re (g\g)^ = p(a',t)2Re (g\g)+ a , (33) 
then p(cr, t) is the stationary solution for (EJ2). 
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Let us investigate Q52p, ( |3"BD for the equilibrium state of the field. In this case 



2Re (g\g)- a , = 2vr J dk \g(k)\ 2 5(u(k) + s a , - e^^—^ 

2tt j dk \g(k)\ 2 5(uj(k) + e a > - e a ) 
2Re (g\g)t a , = ■ 



Crrrr' 



1 — e -P(£<T-s a i) \ — e -/3(e CT -e CT /) 

Crrrr' 



e P(e a -e a i) _ I 

The equation ([32]) takes the form 

U», tV" = E ^^^S^dff^ - *>«*-) (34) 

Let us note that C aa i are non-zero (and therefore positive) only if denominators in ( p4[ ) are positive 
and G a i a are non-zero only if the corresponding denominators are negative. 

If the system possesses decoherence then C aa * , C a i a are non-zero and the solution of equation 



for t — ► oo tends to the stationary solution given by the detailed balance condition (|33 ) 

1 — e -P{z<T-£<?i) eP(e<T-z a i) — 1 



or 



p(a,t)e^ =p{a',t)e^' 

This means that the stationary solution ([33]) of (|32| ) describes the equilibrium state of the system 

p(a, t) 



For a system with decoherence the density matrix will tend, as t — * oo, to the stationary solu- 
tion (|33|) of (|32]). In particular, as t — > oo, the density matrix collapses to the classical Gibbs 
distribution. 

The phenomenon of a collapse of a quantum state into a classical state is connected with the 
quantum measurement procedure. The quantum uncertainty will be concentrated at the degrees 
of freedom of the quantum field and vanishes after the averaging procedure. One can speculate 
that the collapse of the wave function is a property of open quantum systems: we can observe 
the collapse of the wave function of the system averaging over the degrees of freedom of the 
reservoir interacting with the system. Usually the collapse of a wave function is interpreted as a 
projection onto a classical state (the von Neumann interpretation). The picture emerging from 
our considearations is more general: the collapse is a result of the unitary quantum evolution 
and conditional expectation (averaging over the degrees of freedom of quantum field). This is 
a generalization of the projection: it is easy to see that every projection P generates a (non 
identity preserving) conditional expectation Ep(X) = PXP, more generally a set of projections 
Pi generates the conditional expectation 
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but not every conditional expectation could be given in this way. 

We have found the effect of the collapse of density matrix for pit) = (UtpU?), where U t = 
liniA^o e °e~ is the stochastic limit of interacting evolution. The same effect of collapse will 
be valid for the limit of the full evolution e~ ltH , because the full evolution is the composition 
of interacting and free evolution. The free evolution leaves invariant the elements of diagonal 
subalgebra and multiplies the considered above nondiagonal element |cr')(er| by the oscillating 
factor e rf ( £ <7'~ £<T ). Therefore for the full evolution we get the additional oscillating factor, and the 
collapse phenomenon will survive. 



5 Control of coherence 

In this section we generalize the approach of || and investigate different regimes of qualitative 
behavior for the considered model. 

The master equation (|32| ) at first sight looks completely classical. In the present paper we 
derived this equation using quantum arguments. Now we will show that (j32|) in fact describes a 
quantum behavior. To show this we consider the following example. 

Let us rewrite ( |3^) using the particular form (|1^) of ((?|5') ± . Using (fL5|), (|16|) we get 



j t p(a,t)=J2^Jdk\g(k)\^(N(k) + l) 



[p(a',t)5(u(k) +e a - e a ,)\(a\D\a')\ 2 - p(a,t)5(u(k) + e a , - e a )\(a'\D\a)\ 2 ) + 
+N(k)(p(a\t)5(u(k) + e a ,-e a )\{a'\D\o)\ 2 -p(o,t)5(u(k)+ (35) 



The first term (integrated with N(k) + 1) on the RHS of this equation describes the emission 
of bosons and the second term (integrated with N(k)) describes the absorption of bosons. For 
the emission term the part with N(k) describes the induced emission and the part with 1 the 
spontaneous emission of bosons. 

Let us note that the Einstein relation for probabilities of emission and absorption of bosons 
with quantum number k 

probability of emission N(k) + 1 
probability of absorption N(k) 

is satisfied in the stochastic approximation. 

The formula (|35| ) describes a macroscopic quantum effect. To show this let us take the spectrum 
of a system Hamiltonian (the set of system states £ = {&}) as follows: let £ contain two groups Si 
and E 2 of states with the energy gap between these groups (or, for simplicity, two states o\ and <T2 
with Ecr 2 > e ai ). This type of Hamiltonian was considered in different models of quantum optics, 
see for review [[Hj (for the case of two states we get the spin-boson Hamiltonian investigated in 
using the stochastic limit). Let the state (■) of the bosonic field be taken in such a way that 
the density N(k), of quanta of the bosonic field, has support in a set of momentum variables k 
such that 

< uj{k) < ujq < \e ai — £<t 2 |) k G suppA^(Zc) (36) 
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This means that high-energetic bosons are absent. It is natural to consider the state (•) as a sum 
of equilibrium state at temperature and non-equilibrium part. Therefore the density N(k) 
will be non-zero for small k because the equilibrium state satisfies this property. 

Under the considered assumption fl36|) the integral of 5-function S(uj(k) +e ai —s a2 ) with N{k) 
in (|35D equals identically to zero. Therefore the RHS of (j35|) will be equal to 



^2tt / dk\g(k)\ 2 ( p(a' ,1)5(10^) +e a - E a ,)\{a\D\o')\< 

-p(a,t)6(u(k) +e a , - e a )\(a'\D\a) 

It is natural to consider this value (corresponding to the spontaneous emission of bosons by the 
system) as small with respect to the induced emission (for N(k) » 1). In this case the density 
matrix p(a, t) will be almost constant in time. This is an effect of conservation of quantum 
coherence: in the absence of bosons with the energy ui{k) equal to e ai — e U2 the system cannot 
jump between the states o\ and a% (or, at least, this transition is very slow), because in the 
stochastic limit such jump corresponds to quantum white noise that must be on a mass shell. 

At the same time, the transitions between states inside the groups Si and £2 are not forbidden 
by (|36|), because these transitions are connected with the soft bosons (with small k) that are present 
in the equilibrium part of (■). In the above assumptions equation ([35]) describes the transition of 
the system to intermediate equilibrium, where the transitions between groups of states £1 and E 2 
are forbidden. 

If the state (•) does not satisfy the property (|36|), then the system undergoes fast transitions 
between states o\ and 0%. We can switch on such a transition by switching on the bosons with 
the frequency uj{k) = e a , 2 — e ai . 



In conclusion: equation (|35| ) describes a macroscopic quantum effect controlled by the distri- 



bution of bosons N(k) which can be physically controlled for example by filtering. 



6 The Glauber dynamics 

In the present section we apply the master equation (|32"D to the derivation of the quantum extension 
of the classical Glauber dynamics. The Glauber dynamics is a dynamics for a spin lattice with 



nearest neighbor interaction, see |!T[, [12]. We will prove that the Glauber dynamics can be 
considered as a dynamics generated by the master equation of the type ([52) derived from a 
stochastic limit for a quantum spin system interacting with a bosonic quantum field. 

We take the bosonic reservoir space T corresponding to the bosonic equilibrium state at tem- 
perature f3~ l . Thus the reservoir state is Gaussian with mean zero and correlations given by 

(«•(*)«(*')> = ^3i^-*0 

For simplicity we only consider the case of a one dimensional spin lattice, but our considerations 
extend without any change to multi-dimensional spin lattices. 

The spin variables are labeled by integer numbers Z, and, for each finite subset ACZ with 
cardinality |A|, the system Hilbert space is 

T~Ls = T~La = ®reAC 2 
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and the system Hamiltonian has the form 



Hs = H A = — - ^ JrsVr-Vs 
r,s€A 

where ejf, a* are Pauli matrices (r G A) at the r-th site in the tensor product 

<r* = l®---<8)l<g»cr < <g»l®-.-<8»l 

For any r, s G A 

J-rs Jsr £ -R; "rr 

We consider for simplicity the system Hamiltonian that describes the interaction of spin with the 
nearest neighbors (Ising model): 

J rs Jr,r+1 

The interaction Hamiltonian Hi (acting in Tis <8> J~) has the form 

Hi = £ < ® V>G?r), V>(<?) = + A*(</), A(g) = f dkg(k)a(k), 
where ip is a field operator, is a smeared quantum field with cutoff function (form factor) 

g(k). 

The eigenvectors \a) of the system Hamiltonian Hs can be labeled by spin configurations a 
(sequences of ±1), which label the natural basis in Hs consisting of tensor products of eigenvectors 
of o z r (spin up and spin down vectors \e r ), corresponding to eigenvalues e r = ±1) 

|c) = <8VeA|£r) 

In the present section we denote e r the energy of the spin at site r, and denote as E(a) the energy 
of the spin configuration a 

r,seA 

The action of the operator a® on the spin configuration a is defined using the action of 
on the corresponding eigenvector \a): so the operator flips the spin at the r-th site in the 
sequence a (i.e. it maps the vector \e r ) in the tensor product into the vector | — e r )). From the 
form of Hs and Hj it follows that, in (|3"2|), the matrix element (a\D\a') of any two eigenvectors, 
corresponding to the spin configurations er, er', will be non-zero only if the configurations cr, a 1 
differ exactly at one site. If the configurations a, a' differ exactly at one site then (a\D\a') = 1. 

The (classical) Glauber dynamics will be given by the master equation for the density matrix 
laying in the algebra of spectral projections of the system Hamiltonian ( j32[) 

±p{a,t) = J2(p«*,t) (2Re (g\g)-^ a + 2Re {g\g)+ a> ) - 

al reA V 

-p(a,t) (2Re {g\g)%^ + 2Re {<j\g)^ ff )) (37) 
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that gives the Glauber dynamics of a system of spins, see [IT|, [12]. Here 

2Re (g\g)~ a * a = 2n J dk \g{k)\ 2 5{uj{k) - J r _i >r e r _i - Jr.r+iSr+i) 1 _ ( 38 ) 

and analogously all the other (glg)^. 

Up to now we have investigated the dynamics for the diagonal part of the density matrix. 
The master equation for the off-diagonal part of the density matrix ( p5|) will give the quantum 
extension of the Glauber dynamics. We consider now this off-diagonal part: 

^2p(jjl,u, t)\fi)(u\ 



From (25) , (|28) we obtain the equation for the off-diagonal elements of the density matrix 

-j^p(p, v, t) = A^p(p, v, t) (39) 

A i»> = E ( ilm (9\g)»,<T*» - ilm {g\g)v,a*v - * Im (g\g) + ilm (g\g)i*v,v 

-Re (g\g)~p*p - Re {g\g)~ a> - Re (g\g)t^ 4l - Re {g\g)% v ,}j (40) 

Equations (|3~7D , (|3~9|), (|40|) describe the quantum extension of the classical Glauber dynamics (|37|) . 
As it was already noted in section 4, the coefficient in (|40| ) is proportional to |A| (the number 
of particles in the system). Due to the summation on r e A the coefficient will diverge for 
large |A| (the real part of A^ will tend to — oo). Therefore the density matrix will collapse to the 
diagonal subalgebra (the classical distribution function) very quickly. 

Let us consider now the particular case of one dimensional system with translationally invariant 
Hamiltonian: 

The translationally invariant Hamiltonian does not satisfy the generic non degeneracy conditions 
on the system spectrum that we have used in the derivation of equations fl37|) , fl39|) and therefore 
we cannot apply these equations to describe the dynamics for this Hamiltonian. 

However in the translation invariant one-dimensional case we can investigate these equations 
by direct methods. 

In this case the (fl , |fl l ) ± , given by fl38|), are non-zero only if e r _i = e r+ \ = 1, and we get for (|38|) 

2Re (glg)-^ = 2njdk \g(k)\ 2 6(u(k) - 2J) ^ _ L— = . _ ^ (41) 

Therefore for one-dimensional translation invariant Hamiltonians we get for (|3"TD, compare 
with 0,@ 

jA°> *) = T3^w ( E {e-Wptf*, t) - P (a, t)) + 

+ E (p(a r V,t)-e- 2 ^p(a,t))) (42) 

reA;E(cr)<E«cr) ' 
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The detailed balance stationary solution of fl4"2] ) satisfy the following: for two spin configurations 
a, a^a that differ by the flip of spin at site r the energy of corresponding configurations differ by 
2 J. The expectation p(/i), fi = <J, a^a of configuration with the higher energy will be e~ 2/3J times 
less. 

For the off-diagonal part of the density matrix for the case of one-dimensional translation 
invariant Hamiltonian the terms in the imaginary part of ( f4"0"D cancel and using ( f4"T| ) we get for 

^ / 2C 2C \ _ ^ 1 + e~ 2/3J 

- I i _ „-2/3J + P 2/3J _ 1 / ^ 1 _ p-2/3J 

reA \ le e reA le 

This sum, over r, of equal terms diverges with |A| — > oo. Therefore the off-diagonal elements 
of the density matrix that satisfy (|39"D will decay very quickly and for sufficiently large t, |A| the 
dynamics of the system will be given by the classical Glauber dynamics. 

For the master equation considered above we used the master equation for generic (non- 
degenerate) Hamiltonian. This gives us the Glauber dynamics. But the translation invariant 
Hamiltonian is degenerate. Therefore in the translation invariant case we will get some general- 
ization of the Glauber dynamics. To derive this generalization let us consider the general form 
of the master equation. For the considered spin system this gives 

^2—p(fi,u,t)\ii)(u\ = ]TpGu,z/,t) E i l (\v)( p E(v)<P»-tfv\ - \P EM a x a P^a^)(u\) + 
C / 1 

, ^E( M )-jO r »(^M-W^|-o (\^){ P E{ V )<PE{u)-J^\ + \P E {n)P X aP E {»)-jVtlj)(v\ 



-i {\p){P E {u)<y x a Pu + <y x v\ - |Pb ((U )<P m+ (7 6 »(i/|) + 

C 



D 2f3J 



j(\PE(ri+jVaV)(PE(v)+jtfv\ ~ \ (\H)(PE(V)V X PE(U)+J° X V\ + (" 

(43) 

Here C is given by fl4"T|), operator P E u) is a projector onto the states with the energy E(fi), 
operator P u _ is given by 

P„_ = Im (g\g)Z 1 PE(u)-j + Im (g\g) Q P E{u) + Im (g\g)^ P E ( v )+j 

Im (g\g)- = P.P. f dk\g(k)\ 2 —^ ^tttt, a = -1,0,1 

For the operator P, + we get the analogous expression 

P v+ = Im (g\g)tiPE{v)-j + Im {g\g)~% P E (y) + Im Pe(u)+j 
with the coefficients ((? 1^)^: 



The equation (43) gives the quantum generalization of the Glauber dynamics. The matrix elements 
p(/i, v, t) of the density matrix corresponding to the states /x, v with different energies will decay 
quickly. But for the translation invariant Hamiltonian there exist different /i, v with equal energies. 
Corresponding matrix element will decay with the same speed as the diagonal elements of the 
density matrix. Moreover one can expect non-ergodic behavior for this model. Therefore the 
generalization (f43|) of the Glauber dynamics is non-trivial. 
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6.1 Evolution for subalgebra of local operators 

In this section to compare with the results of 0] we consider the dynamics of spin systems, 
described in the previous section, for Hamiltonian with non necessarily finite set of spins A but 
for local observable X. 

The observable X is local if it belongs to the local algebra, that is UHF-algebra (uniformly 
hyperfinite algebra) 

A = (J A A 

A is finite 

where A\ is the *-algebra generated by the elements 

®iX h Xi = 1 for i g A 
Consider now the action of #0 on local X 

MX) =££(-i Im ( 9i \ 9j UX, El (A) E u (Dj)] + 1 Im ( 9i \ 9j )+[X, E w (A) K 

+2Re {g t \g 3 )Z (k (A) XE U (Dj) - ±{X, E* u (A) A (D 3 )} 

+2Re ( 9i \ 9j )+ (E u (A) XE: (Dj) - ±{X, E w (A) K Pi)})) (44) 
For E u (A) we get 

Eu (A) = £ PE{r)-uDiP E { r ) = 
E(r)eF u 

= 1 ® |e f _i) (e 4 _i| ® I -e f ) (e*! ® (e»+i| ® 1 + 1® | <— ® |e 4 )<— e<| ® | -ef+i) (-£ i+ i| ® 1; 

(45) 

with the frequency w of the following form 

Therefore the operator A, given by (|45]) is local and moreover, corresponding map #0 given by 
(|44| ) maps .4 into itself. 

The formula (|45D explains the physical meaning of the operator A (A)- For positive a; this it 
flips the spin at site % along the direction of the mean field of its neighbors (for negative u it flips 
the same spin into the opposite direction). 
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